
*Graph Settings
grstyle clear
set scheme s2color
grstyle init
grstyle set plain, box
grstyle color background white
grstyle set color Set1
grstyle yesno draw_major_hgrid yes
grstyle yesno draw_major_ygrid yes
grstyle color major_grid gs8
grstyle linepattern major_grid dot
grstyle set legend 4, box inside
grstyle color ci_area gs12%50

******Descriptive statistics
***Table 2: Average Treatment Effects

ttest Q134_combined, by(Q134_treatment)

**Figure 1: Responses by Partisanship and Treatment Status
graph drop _all
cibar Q134_combined, over1(Q134_treatment) over2(OppositionPartisan1) barlabel(on) gap(0) graphopts(name(graph_2, replace))
ttest Q134_combined if Partisanship==1, by(Q134_treatment)
ttest Q134_combined if Partisanship==2, by(Q134_treatment)


***Figure 2: Responses by Nonpartisanship
graph drop _all
cibar Q134_combined if Partisanship==0, over1(Q134_treatment) barlabel(on) gap(0) 
ttest Q134_combined if Partisanship==0, by(Q134_treatment)


***Figure 3: Responses by Trust in the President for Nonpartisans
graph drop _all
cibar Q134_combined if Partisanship==0, over1(Q134_treatment) over2(TrustBiyaDichotomous) barlabel(on) gap(0) graphopts(name(graph_1, replace)) 
ttest Q134_combined if Partisanship==0 & TrustBiyaDichotomous==0, by(Q134_treatment)
ttest Q134_combined if Partisanship==0 & TrustBiyaDichotomous==1, by(Q134_treatment)

***Appendix C: Comparing Respondents to Nonrespondents
ttest Gender, by(responded)
ttest AGE, by(responded)
ttest Education, by(responded)
ttest SES_Factor, by(responded)
ttest Q98, by(responded)
ttest Q100, by(responded)
ttest News_Factor, by(responded)
ttest Freedom_Factor, by(responded)

ttest Gender if Partisanship==0, by(responded)
ttest AGE if Partisanship==0, by(responded)
ttest Education if Partisanship==0, by(responded)
ttest SES_Factor if Partisanship==0, by(responded)
ttest Q98 if Partisanship==0, by(responded)
ttest Q100 if Partisanship==0, by(responded)
ttest News_Factor if Partisanship==0, by(responded)
ttest Freedom_Factor if Partisanship==0, by(responded)


***Appendix D: Power Analysis
power twomeans 0.5, diff(0.05 0.1 0.15 0.2 0.25 0.3 0.35) sd(0.49) nratio(1) init(100) graph

*Full Sample
power twomeans 0.5, diff(0.04) sd(0.50) n(100 200 400 800 1600 2400 4800)graph
power twomeans 0.5, diff(0.036) sd(0.49) n(583)
power twomeans 0.5, diff(0.036) sd(0.49) nratio(1)
*Nonpartisans
power twomeans 0.5, diff(0.11) sd(0.50) n(100 200 400 800 1600 2400 4800)graph
power twomeans 0.5, diff(0.112) sd(0.50) n(250)
power twomeans 0.5, diff(0.112) sd(0.50) nratio(1)
*Regime Partisans
power twomeans 0.5, diff(0.09) sd(0.50) n(100 200 400 800 1600 2400 4800)graph
power twomeans 0.5, diff(0.09) sd(0.48) n(174)
power twomeans 0.5, diff(0.09) sd(0.48) nratio(1)
*Opposition Partisans
power twomeans 0.5, diff(0.01) sd(0.50) n(100 200 400 800 1600 2400 4800)graph
power twomeans 0.5, diff(0.01) sd(0.49) n(148)
power twomeans 0.5, diff(0.01) sd(0.49) nratio(1)
*Nonpartisans/High Trust
power twomeans 0.5, diff(0.02) sd(0.48) n(100 200 400 800 1600 2400 4800)graph
power twomeans 0.5, diff(0.02) sd(0.48) n(119)
power twomeans 0.5, diff(0.02) sd(0.48) nratio(1)
*Nonpartisans/Low Trust
power twomeans 0.5, diff(0.24) sd(0.50) n(100 200 400 800 1600 2400 4800)graph
power twomeans 0.5, diff(0.24) sd(0.48) n(122)
power twomeans 0.5, diff(0.24) sd(0.48) nratio(1)


***Appendix E: Randomization Inference P-Values
ritest Q134_treatment _b[Q134_treatment], reps(1000) kdensity: regress Q134_combined Q134_treatment if Partisanship==0
ritest Q134_treatment _b[Q134_treatment], reps(1000) kdensity: regress Q134_combined Q134_treatment if Partisanship==1
ritest Q134_treatment _b[Q134_treatment], reps(1000) kdensity: regress Q134_combined Q134_treatment if Partisanship==2
ritest Q134_treatment _b[Q134_treatment], reps(1000) kdensity: regress Q134_combined Q134_treatment
ritest Q134_treatment _b[Q134_treatment], reps(1000) kdensity: regress Q134_combined Q134_treatment if Partisanship==0 & TrustBiyaDichotomous==1
ritest Q134_treatment _b[Q134_treatment], reps(1000) kdensity: regress Q134_combined Q134_treatment if Partisanship==0 & TrustBiyaDichotomous==0

**Appendix F: Pre-Existing Knowledge of SCT
graph drop _all
cibar Q134_combined, over1(Q134_treatment) over2(Partisanship) over3(heard_SCT) barlabel (on) gap(0)
ttest Q134_combined if heard_SCT==0 & Partisanship==0, by(Q134_treatment) 
ttest Q134_combined if heard_SCT==0 & Partisanship==1, by(Q134_treatment) 
ttest Q134_combined if heard_SCT==0 & Partisanship==2, by(Q134_treatment) 

ttest Q134_combined if heard_SCT==1 & Partisanship==0, by(Q134_treatment) 
ttest Q134_combined if heard_SCT==1 & Partisanship==1, by(Q134_treatment) 
ttest Q134_combined if heard_SCT==1 & Partisanship==2, by(Q134_treatment) 


***Appendix G: Responses by Trust in the President
graph drop _all
cibar Q134_combined, over1(Q134_treatment) over2(TrustBiyaDichotomous) barlabel(on) gap(0) graphopts(name(graph_1, replace)) 
ttest Q134_combined if TrustBiyaDichotomous==0, by(Q134_treatment)
ttest Q134_combined if TrustBiyaDichotomous==1, by(Q134_treatment)


***Appendix I: Alternative Scenario _TBD
graph drop _all
cibar Q132_combined, over1(Q132_treatment) over2(OppositionPartisan1) barlabel (on) gap(0) graphopts(name(terror, replace))
ttest Q132_combined if Partisanship==1, by(Q132_treatment)
ttest Q132_combined if Partisanship==2, by(Q132_treatment)

graph drop _all
cibar Q132_combined if Partisanship==0, over1(Q132_treatment) barlabel(on) gap(0) 
ttest Q132_combined if Partisanship==0, by(Q132_treatment)

graph drop _all
cibar Q132_combined if Partisanship==0, over1(Q132_treatment) over2(TrustBiyaDichotomous) barlabel(on) gap(0) graphopts(name(graph_1, replace)) 
ttest Q132_combined if Partisanship==0 & TrustBiyaDichotomous==0, by(Q132_treatment)
ttest Q132_combined if Partisanship==0 & TrustBiyaDichotomous==1, by(Q132_treatment)
